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ABSTRACT 

Context. The mixing length theory (MUT) used to compute the temperature gradient in superadiabatic layers of stellar (interior 

and atmosphere) models contains in its standard form 4 free parameters. Three parameters are fixed a priori (and define what we 

denote as the MUT 'flavour') whereas one (the so-called mixing length) is calibrated by reproducing observational constraints. The 

'classical'Bohm-Vitense flavour is used in all modern MUT-based stellar model computations and, despite its crude approximations, 

the resulting Teff scale appears - perhaps surprisingly - remarkably realistic, once the mixing length parameter is calibrated with a 

solar model. 

Aims. Model atmosphere computations employ parameter choices different from what is used in stellar interior modelling, raising the 

question of whether a single MLT flavour and mixing length value can be used to compute interiors and atmospheres of stars of all 

types. As a first step towards addressing this issue, we study whether the MLT flavour (the so-called ML2) and mixing length choice 

that have been proven adequate to model white dwarf atmospheres, is able to provide, when used in stellar models, results at least 

comparable to the use of the 'classical'Bohm-Vitense flavour. 

Methods. We have computed solar models and evolutionary tracks for both low- and intermediate-mass Population I and II stars, 

adopting both solar calibrated Bohm-Vitense and ML2 flavours of the MLT in our stellar evolution code, and state-of-the-art input 

physics. 

Results. The two sets of models provide consistent results, with only minor differences. Both calibrations reproduce also the T^ff of 

red giants in a sample of Galactic globular clusters. The ML2 solar model provides a mixing length about half the value of the local 

pressure scale height, thus alleviating - but not eliminating - one of the well known inconsistencies of the MLT employed in stellar 

models. This mixing length is also consistent with the value used in white dwarf model atmosphere computations. 

Key words, convection - globular clusters; general - stars: interiors, evolution - Sun: general - turbulence 



' 1. Introduction 

^^ The equations still largely employed to determine the value of 
;_j the temperature gradient in the superadiabatic regions of stel- 
Ci I lar convective envelopes and atmospheres are based on the so 
called "mixing length theory" (MLT). The MLT in its original 
form is a simple, local , time independen t model, firstly applied 
to stellar modelling bv lBiermannI (|1932|) . Its description of con- 
vective transport considers the motion of 'average' convective 
cells, all with the same physical properties at a given radial lo- 
cation r within the star. These convective cells have a mean size 
I {I = umltHp is the so-called 'mixing length', where umlt 
is a free parameter and Hp is the local pressure scale height 
Hp = P/gp) also equal to their mean free path, and an average 
convective speed Vc at a given value of r. There is in principle 
no reason why umlt should be kept constant when consider- 
ing stars of different masses and/or chemical composition and/or 
at different evolutionary stages; even within the same star umlt 
might in principle vary from layer to layer, although stellar mod- 
els usua lly keep aMir consta nt throughout the convective zone 
(but see lPedersen et al.ll 19901 for some low-mass stellar models 
computed with varying qmlt within the convective envelope). 



Keeping umlt constant still means a change of the cells' mean 
size and free path due to the variation of Hp. 

The value of aj^LT is usually calibrated by reproducing some 
empirical constraint. In case of stellar models, it is typically the 
calibration of a theoretical solar model that determ ines the value 
of a MLT ■ I n all major ste llar model databases (e.g.'Girardi et aL 
200C; J^denBerg et all 120001: lYi et aljllooit [Pietrinferni et ak 
20041 : [Potter et aI.I I2007) this solar calibrated umlt is kept fixed 
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in all evolutionary calculations of stars of different masses and 
chemical compositions. There are 3 additional free parameters 
entering the MLT that are generally fixed a priori, before the 
calibration of omlt- They will be denoted here as a, b, c follow- 
ing the formalism of iTassoul et al.l (Il990l) . These three parame- 
ters plus a MLT enter the equations that determine Vc, the average 
convective flux Fc and the convective efliciency F (defined as 
the ratio between the excess heat content of a raising convective 
bubble just before its dissolution, and the energy radiated during 
its lifetime) at each point within the convective envelope. Once 
fl, b, c are fixed, the temperature gradient at a given location in 
the convective envelope of a stellar model and the resulting Teff, 
depend on the value of umlt- On the other hand, when omlt is 
kept fixed and one or more of the other three parameters are var- 
ied, the temperatur e gradients and T^jj are also affected (e.g. 
lHenvevetal1ll965h . 
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Th e 'classical' formulation of the MLT by iBohm-Vitensd 
(Il958l) that is nowadays almost universally used in theoreti- 
cal computations of stellar evolution models (denoted here as 
the BV58 'flavour' of the MLT) employs a specific set of 
choices for the values for a,b and c. Solar calibrations with 
the present generation of stellar input physics provide a value 
of a MLT typically around ~2.0 (e.g. IVan denBerg et al.' '20001: 
Pietrinferni et al. 2004; Bahcall et alj [2006; Dotter et al 20(M 
Morel & LebretonI |2008|) the precise value (for a given set of 
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input physics) being dependent on the choice for the boundary 
conditions at the model surface. 

The situation for model atmosphere computations is differ- 
ent. The model grids perhaps most used to produce synthetic 
spectra and bolometric corrections for evolutionary phases be- 
fore the White Dwarf (WD) stage are the ATLAS9, MARCS 
and PHOENIX o nes. T he recent ATLAS9 grid described in 
iPietrinferni et al.l ( |2004|) employs gMLr^l-25 , while the lat- 
est MARCS models JGustaifsson etall l2008h adopt a diff'er- 
ent value Q'MLr = 1-5. There are some differences between the 
MLT flavours used in thse two sets of models, and also com- 
pared to t he BV58 formalism used in stellar evolution computa- 
tions (see lCastelli eTai]|1997l:lGustafsson et al.ll2008h . ATLAS9 
model s with a much sm aller aMLT =0-5 have also been com- 
puted dHeiter et al.ll2002h . The recent lBrott & Hauschild^ (l2005h 
PHOENIX grid employs a'MLr=2.0 and (as communicated to us 
by our referee) the BV58 flavour of the MLT. 

Yet another flavour of the MLT is employed in calculations 
of model atm ospheres, spectra and bolometric correc tions for 
WD stars (e.g. lBergeron et al.irr995al:lRohrmannll2001h . It is the 
so-called ML2 ( iTassoul et alj 19901) . with its own specific choice 
of a, b and c, dififerent from the case of BV58. iBergeron et alj 
(Il995bl) have shown how for the case of ZZ Ceti stars, the ML2 
with a MLT -Q -6 provides overall consistency between temper- 
ature estimates from both UV and optical spectrum, observed 
photometry, gravitational redshift mass estimates and trigono- 
metric parallax. The exact form of the MLT em ployed in WD 
interi or models seem to be irrelevant (see, e.g. iFontaine et alj 
1200 lb so that effectively one can employ the ML2 in compu- 
tations of both convective atmospheres and envelopes of WDs. 

This variety of choices regarding the 4 free parameters of the 
MLT should probably not be surprising. As already stated be- 
fore, given the very approximate nature of the MLT, one cannot 
a priori expect the same choice of umlt - and possibly also of 
a, b and c - to be appropriate for all evolutionary phases in both 
convective envelope and atmosphere of stars. It is however in- 
teresting, if only for heuristic purposes, to investigate whether 
a single flavour of the MLT with just one choice (or a small 
range of values) of omlt is adequate to model convective en- 
velopes and atmosphere covering all major evolutionary phases. 
Here we present a first step in this direction, where we have se- 
lected as reference choice the ML2. Our aim is to study whether 
a solar calibration with the ML2 provides a value of umlt close 
to what is employed in WD studies, how this 'new' solar model 
compares with a B V58 solar model, and how evolutionary tracks 
computed with the solar calibrated ML2 compare with the BV58 
calibration and with empirical determinations of T^ff for Red 
Giant Branch (RGB) stars in Galactic glob ular clusters. Based on 
the results bv lGough & WeissI (Il976h and lPedersen et al.l (1 19901) 
we should expect that stellar models computed with different 
choices for a, b and c are largely equivalent once aMur is ap- 
propriately recalibrated (although this has not yet been verified 
along the upper RGB) but we definitely do not know whether the 
same value of qmlt employed in WD model atmosphere studies 
is adequate also for stellar modelling with the ML2. 



Table 1. Values of the free parameters (besides omlt) in the 
BV58 and ML2 flavours of the MLT 



We wish to stress that the aim of this note is not to argue for 
the superiority of the MLT over alternativ e and less crude pre- 
script ions for stellar convection (see, e.g. ICanuto & Mazzitellil 
Il99ll) . rather to investigate whether, within the framework of the 
MLT - that in spite of all its intrinsic limitations is still widely 
used in stellar evolution and model atmosphere calculations - 
and with state-of-the-art physics inputs, it is possible to find a 
single combination of parameters suitable for computing inte- 
rior and atmosphere models of all types of stars with superadia- 
batic convective layers. With 'suitable' we mean that its use pro- 
duces r^//, spectral energy distributions and line profiles consis- 
tent with observations. The precise numerical value of umlt we 
derive for a given MLT flavour, is tied to the specific set of input 
physics adopted in the model computation. We employ what we 
believe are among the best possible physically motivated choices 
for these inputs. Inaccuracies in the current generation of inputs 
for stellar models (and model atmosphere) may of course impact 
the calibration of umlt- 

Section 2 summarizes the MLT formalisms and the choices 
of a, b and c in its ML2 and BV58 flavours. A comparison of 
solar stellar models, low mass stars in the metal poor regime and 
intermediate mass models computed with these two prescrip- 
tions for convection is discussed in Sect. 3, together with the 
comparison with Tg/j- estimates of RGB stars in globular clus- 
ters. Conclusions follow in Sect. 4. 



2. MLT formalisms 

As already mentioned in the Introduction, there are three free 
parameters entering the M LT equations (deno ted here as a, b, c 
following the formalism of lTassoul et al.|[l990b besides the mix- 
ing length. These 3 parameters appear (together with umlt) in 
the following equations for Vc, Fc and F 
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where V is the temperature gradient of a rising (or falling) 
element of matter within the convective region, V is the average 
temperature gradient of all the matter at a given level within the 
convective zone (the quantity needed to solve the stellar struc- 
ture equations) and Q = -(dlnp/dlnr)p. Table[T]summarizes the 
B V58 and ML2 choices for a, b, c. ^From these expressions for 
Vc, Fc and F it is possible to obtain a simple algebraic equation 
wh ose solu tion provides the value of V at a given value of r. 

ICox & Giuli ( 1968) provides a widely used implementation 
of the BV58 flavour of the MLT. Here we show how to transform 
their implementation to describe the ML2. 
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The equation for the effective temperature gradient in 
ICox&Giulil (119681) is 



where ^ is defined as 

V,-V 



(4) 



^ = 



v.-v, 



ad 



Once the value of ( is computed, the knowledge of V, (radiative 
gradient) and V„d (adiabatic gradient) provides immediately the 
actual gradient V. The quantities B and oq entering Eq. |4] are 
obtained through the following relationships: 

V,- V = floA(V- V')^^- 
B=[(AVflo)(V,-V,rf)]'/^^ 

To employ Eq.|4]with the ML2 choices of the constants a, b 
and c (Eqs. [1] |2] |3]l, the quantities oq. A, and B calculated with 
the BV58 have to be transformed as follows: 



ao(ML2) = ao(BV58) - 



A(ML2) = A(BV58) 3 V2 



B(ML2) = B(BV58) | — 
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Fig. 1. Run of the superadiabaticity (V-V,,,/) and the ratio of the 
convective to the total energy flux as a function of the radial 
location in the outer layers of the solar convection zone. Solid 
lines represent the ML2 model, dashed lines the BV58 model. 
The dashed dotted line displays the ratio between the local pres- 
sure scale height and the geometrical distance from the top of 
the convective reaion. 



3. Model comparison 

The canonical way of calibrating umlt (together with the ini- 
tial solar H, He and metal mass fractions X, Y and Z) is to 
compute a theoretical solar model that reproduces - at the so- 
lar age - th e observed solar radius , luminosity and {Z/X)q ra- 
tio (see, e.g. iGough & Weiss|[l976l) . We have calibrated a theo- 
retical solar mod e l (incl uding atomic diffusion) as described in 
iPietrinferni et al.l (l2004h . starting from the pre-Main Sequence 
stage. Our adopted stellar evolution code and state-of-the-art in- 
put physics are fully described in Pietrinferni et al. (2004^ and 
references therein) with the exception of the low temperature 
opacities (T < 10"^ K), for which w e used here the updated 
calculations by Ferguson et al.l ( l2005l) that supersede the older 
[Alexander & FergusonI (Il994h tables. The electron conduction 
opacities ha ve also been update d, by employing the new cal- 
culations by ICassisi et al.l ( |2007|) . although they are not impor- 
tant for the solar calibration. Other relevant physics inputs are 
the equation of state by A. Ir wirQ that reproduces very w ell the 
results of the accurate OPAL (Rogers & Navfonov 2002) equa- 
tion of state, as disc ussed in Pietrinferni et al. (2004); the OPAL 
radiative opacities (llglesias & RogersI Il996|) for T > 10^ K , 
and the NACRE nuclear reaction rates (lAngulo et alj Il999h . 
Atomic diffusion ( without the effect o f radiative levitation) is 
treated according to iThoulet al] d 19941) and the surface bound- 
ary condition s were obtained by inte grating the solar T{t) re- 
lat ionship bv iKrishna Swamvl (Il966h. We employed, ag ain as 
in IPietrinferni et al.1 (|2004 . the iGrev esse & Noels' ('1993') solar 
metal mixture, given the strong conflict (Bahcall et al. 2005) be- 
tween inferences fro m helioseismology an d solar models com- 
puted with the recent lAsplund et al] ( |2005|) results. 



http://freeeos.sourceforge.net/ 



The solar calibration with the BV58 and ML2 provides ini- 
tial abundances 7=0.2755 and Z=0.0201 in both cases. The ac- 
tual radius of the bottom of the convective envelope and He 
mass fraction in the envelope were the same in both calibra- 
tions, namely R^: = 0.7158/?o and 7=0.243. The predicted 
value of Z/X for the present Sun is (Z/X)=0.0244 in both 
cahbrations, in ag r eemen t with {Z/X)q - 0.0245 + 0.005 by 
IGrevesse & Noelsl (Il993b . The calibrated values of omlt are, 
respectively, ao_BV58=2.01 and a-QMLa =0.63. Notice the consis- 
tency of ao,ML2 with the best choice for WD atmosphere mod- 
elling. 

To study the properties of superadiabatic regions calculated 
with the ML2, we compare in the following the run of some 
physical variables along the solar convective envelope, with both 
the BV58 and ML2 models. The difference (V - V^^) in the 
outer layers is displayed in Figure [T] and shows how the peak 
of the superadiabaticity is higher and slightly naiTower with the 
ML2 formulation. This stems from the fact that at the location 
of the (V - Vad) peak, the fraction of the total energy flux carried 
by convection is smaller in the ML2 model (see Fig. [Tj) hence 
a larger superadiabatic gradient is needed to satisfy the energy 
transport equation. The ratio between F^ and the total energy 
flux F is also displayed and, as expected, differences between 
the two models are generally small and obviously restricted to 
the superadiabatic region. An interesting quantity also shown in 
Figure[T]is the ratio between the local pressure scale height and 
the geometrical distance from the top of the convective envelope. 
Around the peak of the superadiabatic region the local value of 
Hp is about 3-3.5 times the distance from the surface. This high- 
lights a well known inconsistency (e.g. [Mazzitelli 1999) when 
using the BV58 calibrated values of omlt, given that with a ~ 
2, the mixing length I in the superadiabatic region is about 6- 
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Fig. 2. Run of T vs P (upper panel) and pvsP (lower panel) in 
the outer layers of our solar models. The ML2 results are dis- 
played as solid lines, the BV58 results as dashed lines. 



7 times longer than the distance z from the surface. Use of the 
ML2 alleviates - but does not eliminate - this conflict, producing 
a local I more comparable with z- 

The internal P - T and P - p stratifications are displayed in 
Fig. |2] Diff'erences in the treatment of the superadiabatic layers 
produce slight modifications in the temperature and density val- 
ues at a given P, that disappear when moving either towards the 
surface or towards the inner adiabatic layers. The largest dif- 
ferences in the stratification appear between log(f)~5.10 and 
log(P)~5.15, close to the region of the superadiabatic peak. 
Here, the temperature in the ML2 model increases more steeply 
towards the interior than in the BV58 case, and the p profile 
is almost flat. Finally, Fig. [3] shows the relative sound speed 

(Vs = J~) diff'erence along the entire structure of the two 

models. The values of Vj are essentially identical except in the 
superadiabatic region (Fig. [TJ as expected. However, the maxi- 
mum difference is only ~ 0.4 % (the ML2 sound speed is larger) 
just below the location of the peak of the superadiabatic gradi- 
ent. Moving inwards the ML2 sound speed becomes smaller by 
0.2 % and then equals the BV58 profile. 

After calibrating the solar model we have computed a 
number of metal poor, low-mass models, from the pre-Main 
Sequence to the tip of the RGB, to compare solar calibrated 
Q^MLT models in a regime diff'erent from the solar envelope and 
test predicted RGB Te/f values against the determinations by 
iFrogel et al.l (Il983l) in a sample of Galactic globular clusters. 
Globular cluster RGB T^ff estimates are one of the key data 
to test the treatment of superadiabatic convection in metal-poor 
low-mass stars (see, e.g. Salaris et al. 2002; Ferraro et al. 20Qa, 
and references therein). 

To this purpose, stellar masses with a lifetime of 12- 
12.5 Gyr have been chosen as representative of objects popu- 
lating the Turn Off and RGB of the observed Colour Magnitude 
Diagrams of Galactic globular clusters. We considered the a- 
en hanced (\afFe]-QA) met al mixture and Y - Z relationship in 
the iPietrinferni et al.l (l2006h a-enhanced models, and computed 




-3.0 -2.0 
log(l-r/Rj 

Fig. 3. Relative sound speed difference between the ML2 and 
B V58 models, as a function of the distance from the centre. The 
bottom of the convective region is marked by a vertical dashed 
line. 



models for [Fe/H]=-2.14, -1.62, -1.01, -0.60, corresponding 
to Z=0.0003, 0.001, 0.004 and 0.01, respectively. Sources for 
the input physics are the same as for the solar calibration. The 
boundary co nditions for these models are still obtained by in- 
tegrating the 'Krishna Swam^ (1 19661) solar T(t). The analysis 
by [VandenBerg et al. (2008) has shown that boundary condi- 
tions from their differen tially corrected model atmospheres (see 
IVandenBerg et al.ll2008L for details) that employ the same MLT 
constants as in the stellar interiors, provide almost the same so- 
lar aMLT and also very simil ar Teff a l ong th e RGB at low metal- 
licities as when the Krishna Swamvl d 19661) T{t) relationship is 
used. 

Figuregldisplays the tracks for [Fe/H]=-2. 14 (M = 0.8 Mq) 
and [Fe/H]=-0.60 (M = 0.92 Mq) from the pre-Main Sequence 
to the tip of the RGB. The ML2 and BV58 tracks at the same 
[Fe/H] are almost identical. Differences in T^ff between the 
two sets of tracks begin to appear at the base of the RGB and 
slowly increase towards the RGB tip, when the superadiabatic 
region is located at lower temperatures and lower densities com- 
pared to previous evolutionary phases. In quantitative terms, they 
amount at most to -40 K for [Fe/H]=-2.14, and -50 K for 
[Fe/H] =-0.60. The luminosity of the RGB bump is the same, 
irrespectively of the MLT flavour, and also the amount of He 
dredged up along the RGB is the same. This is a consequence 
of the agreement between the location (in mass) of the bottom 
boundary of the surface convective regions in both sets of mod- 
els, all the way to the tip of the RGB. Evolutionary timescales 
are unaffected by the change of the MLT flavour 

Figure |5] compares also estimates of the mean Tg// at 
Mi,oi--3.Q for the RGBs of a sample of Galactic globular 
clusters by iFrogel et al.l (Il983h . with the corresponding quan- 
tities from our models. Theoretical bolometric luminosities 
have been transformed into bolometric magnitudes assuming 
A/io/,o=4.75, to have the same zero point of Frogel et al. (1983^ 
bolometric magnitudes. T he cluster [Fe/H] values are on the 
ICarretta & GrattonI (Il997h scale. We assigned an 0.1 dex error 
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Fig. 4. Evolutionary tracks of low-mass models for the labelled 
[Fe/H] values and an a-enhanced metal distribution, from the 
pre-Main Sequence to the tip of the RGB, computed with the so- 
lar calibrated B V58 (dashed lines) and ML2 (solid lines) flavours 
of the MLT (see text for details) 



bar to the cluster [Fe/H] values, and a ±100 K error to the tem- 
peratures (as an estimate of re alistic systematic e rrors in the 
RGB temperature scale, see e.g. lAlonso et al.|[l999l) . 




[Fe/H] 

Fig. 5. Estimates of the average RGB T^j-f at M\,oi - -3 for 
a sample of Galactic globular clusters, compared to the corre- 
sponding values from our models, computed with the solar cal- 
ibrated BV58 (dashed line) and ML2 (solid line) flavours of the 
MLT. The assumed typical error bars are displayed in the lower 
left corner (see text for details). 



Fig. 6. Evolutionary tracks for, respectively, a SM© star with 
solar chemical composition (from the subgiant branch to the 
early Asymptotic Giant Branch phase - upper panel) and a 
1.2Mo star with [Q'/Fe]=0.4 and [Fe/H]=-1.62 (from the pre- 
Main Sequence to the base of the RGB - lower panel). As in 
the previous figures, solid lines display the ML2 models, dashed 
lines the BV58 models. 



Both MLT flavours prod uce T^ff values that appear consis- 
tent with lFrogelet al.l (1 1983b data. We computed the differences 
between the predicted and estimated T^ff for each cluster, by 
interpolating (with a cubic spline) among our results for the in- 
dividual cluster metallicites. We did not find any statistical sig- 
nificant correlation between the Teff differences (observations- 
theory) and the corresponding [Fe/H], for both ML2 and BV58 
models. In case of the ML2 models the average diff'erence 
(observations-theory) is equal to -15 K, with arms of 65 K. The 
same difference for the BV58 models is equal to +32 K, with a 
rms of 65 K. The agreement with observations appear good in 
both cases. The smallest difference is for the ML2 results, but 
both s ets of models provide T^j-f consistent with Frogel et al.l 
(1 19831) if one considers the rms dispersion around this mean val- 
ues and the estimated systematic errors on the RGB T^ff scale. 

We remark here that these models, at variance with the so- 
lar calibration, have been computed without including the effect 
of atomic diffusion, who se efficiency in Galactic globular clus - 
ter stars is still debated dGrattonet alj|200ll: iKorn et al.ll2006l) . 
In any case, the efficiency of diffusion does not affect differen- 
tial comparisons along the Main Sequence, because of the same 
mass extension of convective envelopes in both ML2 and BV58 
models (hence the same rate of depletion of surface He and met- 
als in both cases, and the same effect on Teff)- Also, the compar- 
ison with estimated T^fj values of RGB stars is robust, given that 
the location of the RGB in the H-R diagram is largely insensitive 
to the efficiency of diffusion. 

As a final test, we display in Fig.|6]models for a I.IMq star 
with an a-enhanced, Z=0.001 ([Fe/H]=-1.62) metal composi- 
tion, and a 5Mq star with solar composition, with both solar cal- 
ibrated BV58 and ML2 convection. Both objects display during 
their evolution shallow or vanishing surface convective regions, 
and allow us to compare the two MLT flavours in the regime 
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of very thin convective envelopes. The ML2 and BV58 calibra- 
tions produce again identical tracks along the Main Sequence 
and turn off region of the 1.2Mq star, where convection is effi- 
cient in a very shallow outer layer. Differences in T,,// along the 
pre-Main Sequence and RGB are as for the lower mass objects 
discussed above. Evolutionary tracks with the BV58 and ML2 
for the 5Mq solar composition model are also shown in Fig. |6] 
We have displayed the post-Main Sequence evolutionary phases 
where surface convection appears in the models. Differences in 
Teff are again small, reaching at most 50 K along the bright 
RGB and the early Asymptotic Giant Branch. The amount of He 
dredged up during the post-Main Sequence phases is the same 
with both treatments of surface convection, as a consequence of 
equal depths of the convective regions. 

4. Conclusions 

In the previous section we have shown that solar calibrated 
BV58 and ML2 stellar models computed with the input physics 
and boundary conditions described in Sect. [3] are very similar 
There are some differences in the physical structures of the su- 
peradiabatic layers - as expected - together with minor varia- 
tions in Teff, but the extension of the convective regions are the 
same. Effective temperature differences appear essentially at the 
RGB and slowly increase towards the RGB tip, reaching at most 
40-50 K at the low est temperatures. Thi s latter result does not 
conform exactly to iPedersen et al.l (1 19901) conclusion that mod- 
els computed with different choices for the MLT parameters a, 
b and c should be equivalent once aMur is recalibrated, but the 
differences in the RGB Tejf are not large. 

Both ML2 and BV58 solar calibrated models provide a good 
match to Tejj estimates for bright RGB stars in Galactic glob- 
ular clusters. ML2 models appear to perform marginally better 
when observation-theory differences for each cluster are com- 
pared, but the effect of changing the MLT flavour is very small 
compared to both the rms scatter of the individual values, and 
reaUstic estimates of systematic errors in the RGB T^ff scale. 

In solar models, we have compared in some detail the sound 
speed, T, p stratifications in the convective layers. If the ML2 im- 
proves the agreement between predicted and observed helioseis- 
mological properties (e.g. through the compari son of the solar p 
mode excitation rates, see ISamadi et al.ll2006l) is a matter to be 
established. On the other hand, due to its crudeness, we do not 
expect the MLT to be the adequate tool to model local details 
of superadiabatic convection to the degree of accuracy needed 
to satisfy helioseismological constraints on the structure of the 
solar superadiabatic layers. 

A crucial point is that omlt in the ML2 solar calibration ob- 
tained with our code and the adopted state-of-the-art stellar input 
physics is equal to 0.63, a value consistent with the calibration 
preferred in WD model atmosphere calculations, that deal with 
a very different regime of gas properties. The exact value of the 
solar umlt may surely be subject to some variations, when dif- 
ferent sources for the input physics are employed, but we expect 
that with nowadays best (i.e. most accurate and updated) possi- 
ble choices, values around umlt ~ 0.6 should still be found. 

To verify whether this unification of the MLT treatment can 
be pushed even further, it will be necessary - as a next step - to 
compute and test model atmospheres (plus spectra and bolomet- 
ric corrections) for the pre-WD phases with surface convection, 
using the ML2 formalism with omlt ~0.6. These ML2-based 
model atmospheres, when computed and as a further test of this 
unified MLT treatment, should then be also employed to deter- 
mine the boundary conditions of interior models, for studying 



the consistency with results obtained with a (solar) ^(t) integra- 
tion. 
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